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Abstract 

Wc have performed a search for scalar top quark (stop) pair production in the inclusive electron-muon- 
missing transverse energy final state, using a sample of pp events corresponding to 108.3 pb~^ of data 
collected with the D0 detector at Fermilab. The search is done in the framework of the minimal 
supersymmetric standard model assuming that the sneutrino is the lightest supersymmetric particle. 
For the dominant decays of the lightest stop, i — > 6x1^ and t — » hiv, no evidence for signal is found. 
We derive cross-section limits as a function of stop (i), chargino (Xi^)i and sneutrino {v) masses. 
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Supersymmetry (SUSY) Q provides a theoretically 
attractive and coherent picture of the microscopic world 
that retains the standard model's successful description 
of the observed elementary particles and their interac- 
tions. A major consequence of the realization of SUSY 
in nature would be the existence of additional particles 
(sparticles) , with quantum numbers identical to those of 
the elementary particles of the standard model (SM), but 
with spins differing by a half unit. From experimental 
evidence, the sparticle masses also differ from those of 
their SM partners, i.e., SUSY is a broken symmetry, and 
it is expected that the mass spectrum of the sparticles 
has a different pattern than that of the SM. In particu- 
lar, in several SUSY models, the large mass of the top 
quark (mt) induces a strong mixing between the super- 
symmetric partners of the two chirality states of the top 
quark leading naturally to two physical states, ti and 
t2, of very different mass The lightest stop quark ii 
(called i in this Letter) could therefore be significantly 
lighter than the other squarks rendering it a particularly 
auspicious choice for a direct search. 

The production of a pair of stop quarks (it) at the 
Tevatron can proceed through gluon fusion or quark anni- 
hilation. The cross section for such a process depends 
to a large extent only on the stop mass m^, and is 
known at next-to-leading order (NLO) with a precision 
of ±8% [||. The phenomenology of stop decays depends 
on the assumptions of the SUSY model, and this analysis 
is done in the minimal supersymmetric standard model 
(MSSM) D framework with i?-parity g conservation, 
implying that the lightest SUSY particle (LSP) is sta- 
ble. Searches for stop production have already been per- 
formed at the Tevatron assuming that the lightest neu- 
trahno (x?) is the LSP §. 

In this Letter we also search for light stop {rrif < nit) 
production, but assume that the sneutrino (iy) is the LSP. 
Stop searches have been performed under these assump- 
tions at LEP 2 and by the CDF collaboration at the 
Tevatron ||] yielding a mass limit > 123 GeV for 
the lowest allowed sneutrino mass, m^; « 45 GeV, as 
determined at LEP 1 Although these analyses are 
interpreted in the framework of the MSSM, the results 
are largely model independent, depending mainly on the 
masses of the stop and its decay products. 

In the stop mass range probed by the Tevatron, either 
the 2-body decay via a chargino, i — > 6%^, is kinemat- 
ically allowed and thereby dominant, or the chargino 
mediating the decay is virtual and the dominant decay 
mode is t ^ hlu. The three other 3-body decays medi- 
ated by a chargino, i hvl'^ {-^ hvl+x\), i — > bWxi 
and t bH~^Xi7 with subsequent decays Xi^ are 
disfavored In this Letter, the chargino is taken 

either as virtual with a propagator mass of 140 GeV, or 
its mass is varied between its lowest experimental limit 
(« 103 GeV [|ll|) and the maximum value allowed by 
kinematics. The masses of the sneutrinos of all three 
flavors are taken to be equal, except when the channel 



t brUr is assumed to be dominant. 

The experimental signature for decays of a tt pair con- 
sists of two b quarks, two leptons, and missing transverse 
energy (^t)- The variable represents the measured 
imbalance in transverse energy due to the two escaping 
sneutrinos. The leptons can be e, ^ or r, but r leptons 
are considered only if they decay into evD or fxi'D. We 
place no requirements on the presence of jets and use 
only the efil^x signature since it has less background 
than the eelfix or ^^I^t channels. The resulting event 
sample corresponds to 108.3 pb~^ of data collected by 
the D0 experiment at Fermilab during the Run I of the 
Tevatron. 

A detailed description of the D0 detector and its trig- 
gering system can be found in Ref. The data and 
pre-selection criteria are identical to those used in the 
published tt cross section analysis for the dilepton chan- 
nel which includes the selection of events contain- 
ing one or more isolated electrons with > 15 GeV, 
one or more isolated muons with E!^ > 15 GeV, and 
> 20 GeV. is obtained from the vector sum of 
the transverse energy measured in the calorimeter and in 
the muon spectrometer system. Electrons are required 
to have \r]det\ < 1-1, or 1.5 < \r]det\ < 2.5, where rjdet is 
the pseudorapidity (77) defined with respect to the center 
of the detector. Muons must satisfy |77dct| < 1-7. 

The dominant SM processes that provide the e^IpT sig- 
nature are, in order of decreasing importance: i) multi- 
jet processes (called "QCD" in the following) with one jet 
misidentified as an electron and one true muon originat- 
ing from another jet (muon misidentification has negligi- 
ble effects on our final state); ii) Z — > r+r"— s- e^vvvv] 
iii) WW—i eiJLUV] iv) ti^ e^ivDjj; and v) Drell-Yan 
{DY) — > r+T"— + efivvvv. The QCD background was 
determined from data, following the procedure described 
in Ref. [|^. The other backgrounds were simulated and 
reconstructed using the full D0 analysis chain. 

Simulation of the signal is based on Pythia [ p^ , 
using the CTEQ3M ||l^ parton distribution functions 
(PDFs), and the standard hadronization and fragmen- 
tation functions in Pythia. Comphep ||i^ is used to 
generate the 2 and 3-body decays of the stop. Detec- 
tor simulation is performed using the fast D0 simula- 
tion/reconstruction program, which has been checked 
extensively on a reference sample passed through the 
full D0 analysis chain. The ii samples were simulated 
for stop (sneutrino, chargino) masses varying between 
50 (30, 100) and 150 (90, 170) GeV. 

Distributions in the kinematic quantities {E^, Et^, I^t) 
are shown in Fig. |l](a-c). Also shown (d) are the distri- 
butions for the transverse energy of any associated jets, 
defined by a cone algorithm and having Ei^*' > 15 GeV, 
and two additional kinematic quantities in which the 
signal and background display a different response: (e) 
A^^ = \ipe — ip^\, where ipi is the azimuthal angle of the 
lepton £, and (f) Sf^^ = {rje + Vt^l- Based on simulation 
studies, two additional criteria, 15° < A^** < 165° and 
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FIG. 1. Distributions after pre-selection for the total back- 
ground (open histogram), the sum of the total background 
and the expected stop signal for (mp) = 120 (60) GeV 
(shaded histogram) and the data (points) of (a) the transverse 
energy of the electron, (b) the transverse energy of the muon, 
(c) the missing transverse energy, (d) the transverse energy 
of the jets, (e) the difference in azimuthal angle between the 
two leptons, and (f) the absolute value of the sum in 77 of the 
two leptons. 



Y/^^ < 2.0, were applied to improve the signal to back- 
ground ratio in the final sample. 

The expected cross sections for the background pro- 
cesses, the normalized numbers of events passing the pre- 
selection and those passing the final selection are given 
in Table |, and compared to the expected stop signal for 
nij: (mp) = 120 (60) GeV. The eSiciency for selecting the 
signal varies typically between 1% and 4%. The most 
significant sources of uncertainties on the signal are the 
trigger and lepton identification efficiencies (« 12%), the 
stop pair production cross section (8%), the uncertainty 
due to the PDFs (5%) iQ, the effect of the analysis cri- 
teria (6%) and the luminosity (5.3%), which combine to 
approximately 18%. This uncertainty also includes the 
effect of the variation of the SUSY parameters /igusy (the 
higgs-higgsino mass parameter) and p^ . The sys- 
tematic error for the background is about 10%. This 
error is dominated by the uncertainty on the QCD back- 
ground (7%) and on the cross sections for the background 
processes (10-17%). 

The agreement between the number of observed events 
and the expected background leads us to set cross-section 
limits on stop quark pair production. The 95% confi- 
dence level (C.L.) limits are obtained using a Bayesian 
approach ||20|| that takes statistical and systematic uncer- 



Process 


Cross section 
(Pb) 


Events after 
pre-selection 


Events after 
final selection 


"QCD" 


— 


15.1 ± 1.3 


6.7 ±0.5 


-A- — 


1.70 


5.3 ± 1.0 


1.4 ± 0.3 


WW 


0.69 


4.4 ±0.7 


3.3 ±0.3 


tt 


0.40 


2.7 ±0.5 


2.2 ±0.4 


DY-* r+T- 


0.35 


0.18 ±0.04 


0.04 ± 0.02 


Total backg. 




27.8 ±2.7 


13.7 ± 1.5 


Data 




24 


10 


it 


4.51 


17.3 ±3.1 


13.2 ± 2.3 



TABLE I. Cross sections for the background processes, 
the expected numbers of simulated events passing the 
pre-selection and the final analysis criteria, numbers of events 
selected in the e^l^T data sample and the expected stop signal 
assuming (mp) = 120 (60) GeV. 



tainties into account. Assuming that the stop decays via 
a virtual chargino and ma — 50 GeV, any stop mass 
between 73 and 143 GeV is excluded. The CDF collabo- 
ration has also performed a search in the t —^ Mv chan- 
nel 1^, but based on a different signature: large missing 
transverse energy, at least one lepton, one jet identified 
as a & jet, and at least another jet. The CDF and D0 
results are compared in Fig. ^ 

In the MSSM, when the ratio of the two vacuum expec- 
tation values of the Higgs fields is large (tan/3>10), 
the v-T can be substantially lighter than the or the 
f'^, leadingto an enhancement of the decay width for 
i — > brvr [ pOpl| ]. In this case, the absence of signal pro- 
vides a limit on the cross section in this decay channel, 
as shown in Fig. |^ for mo^ = 50 GeV. 

Again assuming lepton universality, more tt production 
limits are shown in Fig. |^ for different values. For a 
fixed value of , the cross-section limit becomes stronger 
with decreasing sneutrino mass, although the difference 
between limits obtained for different mf) decreases for 
high TOj. For mo up to 85 GeV, and for certain values of 
TOj, these are below the expected MSSM cross sections. 



The resulting exclusion contour in the (m^ 



plane 



is displayed in Fig. ^, and compared to those obtained by 
CDF §, LEP 1, and most recently at LEP 2 |||. The 
present analysis places limits at significantly higher m^ 
compared to these results. This is mainly because of the 
higher center of mass energy of the Tevatron compared to 
LEP, and of the choice of a more sensitive signature com- 
pared to CDF. For mo = 45 GeV, the excluded region 
extends up to a scalar top mass of 144 GeV, to be com- 
pared to approximately 123 (98) GeV for CDF (LEP 2). 

The 2-body decay into a b quark and a real chargino, 
i ^ byt , was simulated for m~+ between 100 and 

140 GeV, and the Xi was assumed to decay only into 
Iv, leading to the same final state as i ^ blv. Figure ^ 
shows exclusion contours as a function of m^ and m^+ , 
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■ D0 95% C.L. limit (bxv^j 
A CDF 95% C.L. limit (blvj 
• D0 95% C.L. limit (blv; 




mill <\>i< 2mj 



100 110 120 130 

m \ (GeV) 

FIG. 2. Cross-section limit as a function of mf for 
mc — 50 GeV. The t — > MD results of this analysis are com- 
pared to those of CDF and to the expected NLO cross section 
for three different choices of factorization scale fj,. The renor- 
malization scale is taken to be equal to fi. Also shown is the 
limit obtained in the t hrVr channel for = 50 GeV. 



BR(t^blv)= 100% 




70 90 110 130 

mi (GeV) 



FIG. 4. Excluded regions in the {mi,mp) plane for the 
t blv decay channel in the MSSM. The results of this anal- 
ysis (labelled D0 108 pb~^) are compared to the exclusion 
limits obtained in the t blv decay channel at the Tevatron 
(CDF), and at LEP 2. Also shown is the sneutrino mass limit 
obtained at LEP 1. 
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FIG. 3. Limits on the stop pair production cross section 
as a function of mj, for = 60, 70, 80 and 90 GeV. These 
limits are compared to the expected NLO cross section for 
three different choices of factorization scale /i. 



BR(f^ b Xi) = 100% 

BR(jc|^ 1 V) = 100% 

mr. = 45 GeV 
m~ = 60 GeV 
= 75 GeV 




FIG. 5. Excluded regions in the (mj,m-+) plane for the 

^1 

> hx^ decay channel in the MSSM, for mo = 45, 60 and 
75 GeV. These results are compared to the exclusion limit 
obtained at LEP 2. 
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assuming mo = 45, 60 or 75 GeV. They are compared to 
the exclusion Hmit obtained at LEP 2 assuming unifica- 
tion of the gaugino masses and decay of the chargino via 
aH^*0. 

In conclusion, our analysis that assumes the v to be 
the LSP places new limits on the stop mass. Assuming 
lepton universality and a virtual intermediary chargino, 
the excluded region at 95% C.L. extends up to a scalar 
top mass of 144 (130) GeV for mo = 45 (85) GeV. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the Depart- 
ment of Energy and National Science Foundation (USA) , 
Commissariat a L'Energie Atomique and CNRS/Institut 
National de Physique Nucleaire et de Physique des Par- 
ticules (France) , Ministry for Science and Technology and 
Ministry for Atomic Energy (Russia), CAPES and CNPq 
(Brazil), Departments of Atomic Energy and Science and 
Education (India), Colciencias (Colombia), CONACyT 
(Mexico), Ministry of Education and KOSEF (Korea), 
CONICET and UBACyT (Argentina), The Founda- 
tion for Fundamental Research on Matter (The Nether- 
lands), PPARC (United Kingdom), Ministry of Educa- 
tion (Czech Republic), and the A. P. Sloan Foundation. 



C 15, 1 (2000 ) and 2001 off-year pa rtial update for the 



2002 edition), ittp : //pdg. Ibl . gov/ 



[10] W. Porod, Phys. Rev. D 59, 095009 (1999). 
[11] LEP SUSY Working Group, ALEPH, DELPHI, L3 and 

OPAL Collaborations, LEPSUSYWG/0 1-03.1 (2001), 

http : //lepsusy . web . cern . ch/lepsusy/ . 
[12] D0 Collaboration, S. Abachi et al, Nucl. Instr. and 

Methods A338, 185 (1994). 
[13] D0 Collaboration, S. Abachi et ah, Phys. Rev. Lett. 79, 

1203 (1997). 

[14] D0 Collaboration, S. Abachi et ah, Phys. Rev. D 52, 
4877 (1995). 

[15] T. Sjostrand, Comp. Phys. Commun. 89, 74 (1994); S. 

Mrenna, Comp. Phys. Commun. 101, 232 (1997). 
[16] R. Brock et al, R ev. Mod. Phvs. 6 7, 157 (1995). 
[17] A. Pukhov et al, tiep-ph/9908288 . 
[18] W. Beenakker et al., Nucl. Phys. B515, 3 (1998) 
[19] B. Olivier, Ph. D. Thesis, University of Paris VI (2001) 

(unpublished) . 
[20] I. Bertram et al, Fermilab-TM-2104 (2000). 
[21] W. Porod and T. Wohrmann, Phys. Rev. D 55, 2907 

(1997); A. Datta, M. Guchait and K. K. Jeong, Int. J. 

Mod. Phys. A 14, 2239 (1999); A. Djouadi and Y. Mam- 

brini, Phys. Rev. D 63, 115005 (2001). 
[22] LEP SUSY Working Group, ALEPH, DELPHI, L3 and 

OPAL Collaborations, LEPSUSYWG/0 1-02.1 (2001), 

http : / /lepsusy . web . cern . ch/lepsusy/ . 



* Visitor from University of Zurich, Zurich, Switzerland. 

[1] Y. Golfand and E. Likthman, JETP Lett. 13, 323 (1971) 
D. Volkov and V. Akulov, Phys. Lett. B 46, 109 (1973) 
J. Wess and B. Zumino, Nucl. Phys. B 70, 31 (1974) 
ihid. 78, 1 (1974). 

[2] J. Ellis and S. Rudaz, Phys. Lett. B 128, 248 (1983); M 
Drees and K. Hikasa, Phys. Lett. B 252, 127 (1990) 



[3] W. Beenakker, R. Hopker, M. Spira, |hep-ph/9611232 , 
[4] H. P. Nilles, Phys. Rep. 110, 1 (1984); H. E. Haber and 

G. L. Kane, Phys. Rep. 117, 75 (1985). 
[5] P. Fayet, Phys. Lett. B 69, 489 (1977); G. R. Farrar and 

P. Fayet, Phys. Lett. B 76, 575 (1978). 
[6] D0 Collaboration, S. Abachi et al, Phys. Rev. Lett. 76, 

2222 (1996); D0 Collab., S. Abachi et al, Phys. Rev. 

D 57 589 (1998); CDF Collab., T. Affolder et al, Phys. 

Rev. Lett. 84, 5704 (2000); ihid. 84, 5273 (2000). 
[7] ALEPH Collaboration, R. Barate et al, Phys. Lett. B 

469, 303 (1999); DELPHI Collab., P. Abreu et al, Phys. 

Lett. B 496, 59 (2000); L3 Collab., M. Acciarri et al, 

Phys. Lett. B 471, 308 (1999); OPAL Collab., G. Abbi- 

endi et al, Phys. Lett. B 456, 95 (1999). 
[8] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 

84, 5273 (2000). 
[9] Particle Data Group, D. E. Groom et al, Eur. Phys. J. 



7 



